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TURBULENCE MODELING WORKSHOP 

Edited By R. Rubinstein 1 , C.L. Rumsey 2 * , M.D. Salas’, and J.L. Thomas 4 

Executive Summary. Advances in turbulence modeling are needed in order to calculate high Reynolds 
number flows near the onset of separation and beyond. To this end, the participants in this workshop made the 
following recommendations. (1) A national/intemational database and standards for turbulence modeling assessment 
should be established. Existing experimental data sets should be reviewed and categorized. Advantage should be 
taken of other efforts already underway, such as that of the European Research Community on Flow, Turbulence, 
and Combustion (ERCOFTAC) consortium. Carefully selected “unit” experiments will be needed, as well as 
advances in instrumentation, to fill the gaps in existing datasets. A high priority should be given to document 
existing turbulence model capabilities in a standard form, including numerical implementation issues such as grid 
quality and resolution. (2) NASA should support long-term research on Algebraic Stress Models and Reynolds 
Stress Models. The emphasis should be placed on improving the length-scale equation, since it is the least 
understood and is a key component of two-equation and higher models. Second priority should be given to the 
development of improved near-wall models. Direct Numerical Simulations (DNS) and Large Eddy Simulations 
(LES) would provide valuable guidance in developing and validating new Reynolds-averaged Navier-Stokes 
(RANS) models. Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent 
viable approaches for analysis on a limited basis. Therefore, although computer limitations require the use of RANS 
methods for realistic configurations at high Reynolds number in the foreseeable future, a balanced effort in 
turbulence modeling development, validation, and implementation should include these approaches as well. 

Key words, turbulence, high Reynolds number, separation, modeling, DNS, LES 

Subject classification. Fluid Mechanics 

1. Introduction. In 1996 John Lumley summarized the importance of turbulent flows and our ability to 
calculate them as follows. “Rational design of aircraft, automobiles, nuclear reactors and all sorts of industrial 
mixing and forming process, ... are dependent on an ability to calculate the effects of turbulent transport reliably. 
Unfortunately, we cannot do that. One hundred years of intense effort have brought us very good qualitative 
understanding of turbulent flows in nearly all practical respects, but have not brought us the ability to calculate 
reliably.” [1] While this view by Lumley is correct, there has been considerable progress in the calculation of 
turbulent flows for some restricted classes of flows. The December 1999 report of the Airframe Systems 
Subcommittee of NASA’s Aerospace Technology Advisory Committee (ASTAC) concluded that while great 
success has been achieved with computational fluid dynamics in accurately predicting attached flow, current 
turbulent modeling capabilities are unable to reliably predict separation onset. Indeed, the report considered 
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separation onset as presenting a greater challenge than dealing with massively separated flows. This is perhaps a 
debatable conclusion (see Bradshaw’s presentation in Appendix C). In order to assess our current capabilities and 
future needs for accurate computations of high Reynolds number turbulent separated flows at flight conditions, a 
turbulence-modeling workshop was held in Reno, Nevada on January 12-13, 2001 . The workshop was sponsored by 
NASA Langley Research Center (LaRC) and organized by 1CASE. Approximately 40 technical experts, covering a 
wide range of knowledge, were invited to participate. Ajay Kumar, representing NASA LaRC, opened the workshop 
by establishing its purpose and the expectations he had from the workshop participants. Mark Anderson, Chair of the 
Airframe Systems Subcommittee presented the committee views on current capabilities and shortcomings of 
turbulence modeling. Later that morning and early afternoon, five summary talks were presented. Philippe Spalart, 
Boeing, Brian Smith, Lockheed Martin Aeronautics, and Thomas Gatski, NASA LaRC, presented their own 
perspectives of the state-of-the-art in turbulence modeling, emphasizing high Reynolds number separated flows. 5 
Katepalli Sreenivasan, Yale University, gave an overview of the physics of this flow regime and Roger Simpson, 
Virginia Tech, spoke about issues associated with experimental methods. The rest of the workshop was planned 

around group discussions by the attendees. In order to provide some structure to the discussions, the following three 
topics were chosen: 

a) turbulence modeling for vortical flows, 

b) turbulence modeling for time dependent separated flows, and 

c) turbulence modeling for juncture and mixing flows. 

However, the three topics were not intended in any sense to limit the discussion. The 40 participants were 
divided into three groups and each group was asked to discuss the adequacy of current turbulence models, 
experimental difficulties, numerical issues, and alternative approaches as they related to the three topics above. This 
document summarizes the results of the workshop. 

The organization of this document is as follows. In Sections 2 through 4, summaries of the discussions 
held in each of the three groups are given. These summaries are broken into the following subsections: importance 
of the modeling of turbulence, assessment of current methods, directions for improvement of turbulence models, and 
conclusions and recommendations. Section 5 gives final overall conclusions from the workshop, including an 
assessment of current methods and recommendations for future development. The appendices include an agenda 
from the workshop, a list of participants, and a copy of the workshop presentations. 


5 A fourth talk by Peter Bradshaw was planned, however Bradshaw was unable to attend the workshop. His slides are included in Appendix C. 
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2. Group I - Summary Findings and Recommendations. Facilitators: C.L Rumsey (LaRC) and J.B. 
Anders (LaRC) 

2.1. Importance of the Modeling of Turbulence 

2.1.1. Vortical flows. The following list gives examples of different types of vortical flows of interest to 
the aeronautics/aerospace community: 

• Wing tip vortex 

o Interaction with tail 

o Far downstream 

• Chine vortex 

o Interaction with wing boundary layer, including pressure gradient effects 

• Strake 

• Vortex bursting 

• Fuselage at high alpha and ogive cylinder 

• Vortex generators in boundary layer 

• Internal vortices (separation) 

o Including vortex breakdown 

• Vortex instabilities 

• Flap/junction vortex 

o Mixing enhancers (Chevrons) on engines 

Most vortical flow types on this list fall into one of the following categories: free shear flow or 
vortex/boundary layer interaction. Free shear flow is generally easier to compute, but also tends to be less important 
from the point of view of the aerospace industry (in other words, unless a free vortex comes near a surface, it is not 
so important to compute it accurately for aircraft design). For example, computing wing tip vortices accurately can 
be important when they impinge upon or come near a downstream body (horizontal tail, following aircraft). The 
chine is a protuberance on the outside of the engine nacelle that generates a vortex that can interact with the wing 
boundary layer. The strake (sometimes called strakelet) is a leading edge piece near the wing-body intersection on 
some fighter aircraft, for example, that creates a vortex that passes back over the body and can interact with the 
wing, body, or tail boundary layer. Strake vortices on fighters have been associated with vertical tail buffet (vortices 
that burst induce unsteady loads on the vertical tails as they pass near). 
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It can be more instructive to redefine the above list in terms of physical categorizations, as follows: 

• Free vortex zero pressure gradient (ZPG) 

• Free vortex with pressure gradient 

• Free vortex with and without axial flow 

• Vortex interaction with boundary layer (BL), with and without separation 

• Vortex interaction with shock 

• Interacting vortices (co-rotating and counter-rotating) 

• Smooth body cross-flow separation 

This physical categorization gives a broader representation of the types of vortical flows that can occur. 
Predicting vortex details is not always important from the point of view of the aerospace industry. It depends on the 
case, and tends to be more important when there is an interaction of a vortex with a downstream surface. An 
example was given from Boeing for which loads and moments were accurately predicted in spite of the Spalart- 
Allmaras (SA) turbulence model (which adds too much eddy viscosity in the vicinity of vortices) diffusing the 
vortex prior to the vortex interaction with the tail. Also, it is important to note that engineers in the aircraft industry 
are often looking for accurate predictions of trends , and not absolute levels. On the other hand, the prediction of 
absolute levels of drag is critical to airplane manufacturing and represents one of the most difficult challenges in the 
aerospace industry. 

Generally, Reynolds number is not too important in turbulent free shear flows. However, for free vortex 
flows, the Reynolds number of the vortex-forming device affects the initial vortex formation. If there is turbulence 
decay in a vortex, then viscous transport becomes more important. Far downstream, as the flow becomes quasi- 
laminar, the Reynolds number becomes important. 

2.1.2. Separated/time-accurate flows. The category of separated/time-accurate flows can be broken 
down into several physical categorizations: 

• Curvature (response to normal straining) and pressure gradient 

• Unsteady (hysteresis, time lag) 

• Post-separation physics 

• Post-curvature physics 

• 2-D smooth separation 

• 3-D smooth surface separation 

• 2-D shock-induced separation 

• 3-D shock- induced separation 

• Vortex/BL interaction 

Most separated flows can be categorized in terms of one or more of these physical categories. In addition 
to pressure gradient effects, curvature can affect separation location from a smooth body. Convex surface curvature 
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reduces turbulence whereas concave curvature enhances turbulence. Beyond a separated region, a flow may reattach 
and “recover” from separation. This post-separation physics is often important to compute accurately. For example, 
it is often possible to predict shock location on a wing correctly, but if there is shock-induced separation, its extent 
and possible reattachment location downstream are often not predicted well. Post-curvature physics refers to 
recovery downstream of curvature. Smooth surface separation and shock-induced separation are of particular 
interest to the aerospace community, because of the frequent occurrence of these flow types for aerospace vehicles. 

During this discussion, the issue of Reynolds number scaling was brought up. The comment was made that 
there is no evidence to suggest that Reynolds number scaling is a turbulence modeling issue. 

2.1.3. Juncture and mixing flows. Other types of flows (for example, jets and mixing layers) are 
probably far more important than that of juncture flows from the point of view of the aerospace industry. Accurate 
computation of secondary vortices is probably not important for most typical industrial needs. A prioritized list of 
flows of importance for turbulence modeling is presented in the Conclusions and Recommendations section below. 

2.2. Assessment of Current Methods 

2.2.1. Vortical flows. Often, the problem with capturing vortex interaction effects has less to do with 
the turbulence model than it does with lack of grid refinement in the region where the vortex exists. In this regard, 
automatic grid refinement or adaptation (putting enough grid in the right place) might go a long way toward helping 
to achieve more accurate vortical flow computations. 

Nonetheless, many existing turbulence models are deficient for turbulent vortical flows. To get the details 
right, a turbulence model needs to correctly represent the relationship between stress and curvature. Many models, 
particularly eddy viscosity models (EVM), cannot do this well. Many models erroneously produce eddy viscosity in 
the vicinity of free vortices, which causes the vortices to be excessively diffused above and beyond the effects of 
insufficient grid resolution. 

A general “Pros and Cons” list for turbulence model types is given here, as related to vortical flows (“- 
(minus)” indicates con, and “+ (plus)” indicates pro): 

• EVM - too diffusive 

• EVM+suppression + suppresses diffusion 

- not necessarily at the correct rate (v’w* radial vs. u V axial) 

• ASM/EASM - needs curvature correction 

+ represents normal stress differences 

• RSM - cost/robustness 

+ can reproduce correct behavior 

• LES/DES - cost 

-4- should correctly predict any free shear flow 
(no consensus on this “pro” statement) 
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Key: EVM=eddy viscosity model, EVM+suppression indicates eddy viscosity models with one of many 
available simple fixes that make the model “turn off’ within vortices, ASM=algebraic stress model, 
EASM=explicit algebraic stress model, RSM=Reynolds stress model, LES=large eddy simulation, 
DES=detached eddy simulation. 

Based on the breakdown of vortical flows into physical categorizations given earlier, a table is presented 
here to list the model types that are capable of solving each category. Note that this table should be viewed as a 
framework only. The group did not have the time or all the information necessary to adequately complete it. A 
question mark was used when there was some uncertainty. For example, if a model type has been validated only for 
a single specific case (and some uncertainty as to the model’s validity remains) then that model was assigned “Y?” 
Similarly, limited success in a validation earned a “Y?” and a belief that a model should be capable in spite of its not 
yet being validated also earned a Y? The two items with “N / Y?” are labeled as such because of differing 
opinions given by members of the group. 

Note that DES should work for all cases except boundary layer interaction. The other flows are all free 
shear flows, for which DES defaults to LES. And for any free shear flow, LES should yield good results on a 
sufficiently fine grid, regardless of the Reynolds number (assuming that the spectral content of any inflow boundary 
condition is known, or that the spectral content is not important). 


Table 2.1 



EVM 

EVM with 
suppression 

ASM/EASM 

RSM 

A. Free vortex with zero pressure gradient 

N 

Y 

Y 

Y 

B, Free vortex with pressure gradient 

N 

Y? 

7 

Y 

C. Free vortex with and without axial flow 

N 

7 

? 

Y? 

D. Vortex interaction with BL (with and without 
separation) 

N/Y? 

Y? 

Y with tweak 

Y 

E. Vortex interaction with shock 

N 

7 

7 

Y? 

F. Interacting vortices (co-rotating and counter- 
rotating) 

N 

? 

7 

Y? 

G. Smooth body cross-flow separation 

N/Y? 

Y? 

7 

Y? 


The fixes used in EVM+suppression models (tuning to be sensitized to a curvature parameter such as 
Richardson number) may work well for vortices in which straining in the stream-wise direction is relatively weak. 
However, they do not work in general for vortices interacting with a boundary layer. Also, in more complex 
environments (such as in adverse pressure gradients, in the presence of additional shear, when there is significant 
stream-wise strain, or when there are circulation changes) these simple fixes may not work either. An example of 
one such complex flow is an internal flow in which the shear stresses decay rapidly, but the normal stresses do not, 
so a nearly +ntropic turbulence results downstream. In this region, the turbulent kinetic energy is fairly high. 
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“EVM+suppressiorT may get the representation of the shear stresses right, but it will not be able to compute the 
high turbulent kinetic energy. 

As seen in Table 2.1, RSM is the lowest order model that is generally capable of solving all of the above 
seven flow categories, even though some modeling of terms is required. If a model provides the correct relationship 
between stresses and strains (as implied in the exact Reynolds stress transport equation), then that model should be 
able to yield a good representation of vortices. One pays a price each time a simplification to the equations is made. 
For example, in simplifying from RSM to ASM/EASM, assumptions are made regarding the diffusion term and the 
convective terms for the stresses, which can lead to misrepresentation in certain circumstances. The diffusion term 
can be important in the region of the centerline of the vortex. However, ASM/EASM does capture the interaction 
between stresses and strains as embodied in the stress generation terms. 

In general, there has not been enough CFD validation on many of the above seven flow categories. 
However, validation requires either experimental or full simulation (LES/DNS) data. There are many experimental 
studies that include vortical flows, but the group participants knew of only two direct numerical simulation (DNS) 
studies in this area (both DNS studies are unpublished Ph.D. theses). Clearly, this is an area that could use some 
additional attention. Many of the existing experimental data are very old (for example, Langley facility data). Also, 
many experimental datasets have been taken for delta wings, including vortex breakdown. A thorough survey of 
existing data and its relevance and quality would be helpful. 

2.2.2. Separated/time-accurate flows. Using the physical categorizations for separated/time-accurate 
flows given earlier, the following table lists the model types that are capable of solving each category. As with the 
table above, this table should be viewed as a framework only. Question marks indicate uncertainty, and boxes with 
both N and Y indicate differing opinions offered by the group. 

Note that post-separation physics and post-curvature physics stand out as a challenge to most models, 
including RSM. For shock-induced separation, models often can get the shock location accurately, but the separated 
region itself (downstream of the shock) may be poorly predicted. This behavior is case-dependent. It may in part be 
related to post-separation physics. 

In attempting to assess the capability of existing turbulence models for separated/time-accurate flows, it is 
often difficult to separate whether poor predictions are due to turbulence modeling or other issues. For example, 
transition is often a big problem; if the transition region is unknown in the experiment, then comparing using fully 
turbulent CFD can cause discrepancies. Aeroelasticity, if not accounted for, can also lead to discrepancies. 

Often a given experiment is too complex to be helpful toward isolating specific failings of turbulence 
models. For example, the trapezoidal wing experiment will likely be useful for validation, but is not simple enough 
for improving turbulence models. Simple unit problems are most useful for isolating specific failings of turbulence 
models and guiding modelers toward improvements. (Unit problems isolate a specific aspect of turbulence and are 
also usually geometrically simple so they remove geometric fidelity considerations from the CFD modeling, and are 
easier to grid-converge.) 
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Table 2.2 



EVM 

EVM + 
fixes 

ASM/EASM 

RSM 

A. Curvature (response to normal straining) 
including pressure gradient 

N 

N/Y 

Y 

Y 

B. Unsteady (hysteresis, time lag) 

Y 

Y 

Y 

Y 

C. Post-separation physics 

N 

N 

N 

N/Y? 

D. Post-curvature physics 

N 

N 

N 

N/Y 

E. 2-D smooth surface separation 

N/Y 

Hy 

Y 

Y 

F. 3-D smooth surface separation 

N 

N/Y 

Y 

Y 

G. 2-D shock-induced separation 

Y? 

Y? 

Y 

Y 

H. 3-D shock-induced separation 

7 

? 

7 

? 

I. Vortex-BL interaction 

N/Y? 

Y? 

Y 

Y 


In the area of unit problem experiments, further axisymmetric bump experiments are the type that may be 
helpful for exploring shock-induced separation. However, turbulence models have been very successful in the past 
on simple problems like this. 3-D flows are more challenging. Therefore, in spite of the difficulty inherent in 
defining and carrying out good unit problems for 3-D separation, the turbulence modeling community is definitely in 
need of more “3-D unit problem” experiments. 

2.2.3. Juncture and mixing flows. There is a lot of evidence that full Reynolds stress models (RSM) 
are required for many of these types of flows (for example, horseshoe vortices). It is also well known that nonlinear 
terms are required in a turbulence model for it to be able to compute secondary motions induced by turbulent normal 
stress differences. In other words, if a turbulence model is a linear eddy viscosity model (LEVM), then the turbulent 
stresses are proportional to the strain (Boussinesq assumption) and the model cannot predict turbulent normal stress 
differences. 

2.3. Directions for Improvement of Turbulence Models. Due to time constraints, the group only 
addressed directions for improvement of separated/time-accurate flows. In order to improve turbulence models for 
separated/time-accurate flows, both experiments and turbulence modeling itself must move forward hand-in-hand. 
Turbulence modelers generally like to have, from an experiment: three components of velocity (both mean and 
fluctuating) profiles, temperature profiles, skin friction coefficient, and pressure coefficient. Also, although difficult 
to obtain, some measure of the length scale would be extremely helpful, because the modeling of the length scale is 
currently one of the biggest uncertainties in turbulence modeling. Experimental techniques that are currently 
helpful, and may benefit from further exploitation are: particle image velocimetry (PIV), laser Doppler velocimetry 
(LDV), oil film, liquid crystal, and any other non-intrusive technique. Where experiments are lacking (such as in 
obtaining length-scale information), full simulations may be the only way to move forward. 
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Some goals for improvement in prediction of separated flows are listed here: 

• Increase generality of turbulence model formulation 

• Need naturally good behavior near walls 

• Must include effects of curvature 

• Separation control into modeling 

• Continued exploration of DES 

• Overall, should increase the role of DNS/LES for prediction modeling 

The issue of turbulence modeling implementation is also important. Often, different individuals implement 
a given model differently (for example, different numerical methods, limiters, constants, and/or damping functions 
may be used). Or sometimes, different versions of a given model exist in the literature, yet, when implemented, they 
are referred to by the same name. It is not a trivial task to ensure uniformity, however. In an effort at NASA 
Langley, it took one month for three individuals to modify three existing codes to have identical implementation of 
an EASM model [2]. It is difficult to validate/improve models when such differences exist. (This problem, 
although a much greater problem for more complex models like EASM and RSM, even exists for simple models like 
SA. For example, some major codes in use in industry today employ an unpublished modification to SA that can 
delay the location of transition compared to the published version. Most users are not aware that this modification 
has been employed.) 

2.4. Conclusions and Recommendations. The conclusions and recommendations given here are fairly 
general, not specifically geared toward any one of the flow categories. First of all, it is important to attempt to 
prioritize various flow categories, for which turbulence modeling efforts should be focused. From the point of view 
of the participants, this prioritization is (starting with the most important): 

• Separation (including incipient separation) 

• Vortex flows 

• Jets and mixing layers 

• Unsteady flows 

• Other (Juncture flows. Heat transfer (scalar transport), Flow-induced noise. Compressibility, Cavity 
flows) 

Naturally, this prioritization is subjective. What areas are considered important depend on who is doing the 
prioritizing. However, the top four items in this list seem to the group to represent flows of interest and of 
importance to a great number of people in many disciplines. 

Some specific recommendations follow. First, an effort should be undertaken to make numerical 
implementations for RSM (and other models) efficient (so people will want to use them). There should be standards 
set up to guarantee completeness in reporting details on how turbulence models are implemented (for repeatability). 
Furthermore, details on proper use of models should be published. For example, the guidelines document published 
by ERCOFTAC are a step in this direction [3]. 
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An organized validation effort should be undertaken on a set of simple standard cases. In addition, the 
turbulence modeling community also needs to build up to a more complex set of standard cases. 

After two days of discussion, a lot of attention was given to RSM, as an unspoken goal toward which 
turbulence modeling efforts should be directed in the future. However, the issue was brought up as to whether this 
goal is appropriate or desirable. Is full Reynolds stress modeling an “abyss”? I.e., is it do-able, or is it too much? Is 
it necessary? RSM has been traditionally less robust than simpler models. Can this state of affairs be improved? A 
thorough validation effort is necessary for the existing simpler models before they are discounted. What are specific 
documentable failings that are unambiguously due to turbulence model and not some other factor like grid resolution 
or poor geometric fidelity? (This applies to all turbulence models.) 

The issue of length-scale modeling stood out as one that really needs a lot more focus in the future. This 
will require considerable help from experiments and/or simulations. 

A large concerted effort is needed to evaluate and select from existing experimental databases (such as 
ERCOFTAC). The group did not want to advocate a slew of new experiments when so many old ones exist that 
might serve perfectly well. After a thorough evaluation, some areas may be evident where new or updated 
experiments may be required. Members of the group mentioned some existing experimental databases: trapezoidal 
wing, ROCK wing, and plane swept bump in channel. As a part of compiling an experimental database, a good set 
of unit problem experiments for validation/model improvement is particularly needed. Some suggestions included 
some sort of axisymmetric bump or a modified axisymmetric bump. 

The group made some recommendations for experimental approaches that would be most helpful for 
turbulence modeling validation. These included use of PIV for spatial correlations, LDV, and other non-intrusive 
methods. It is imperative that well-defined boundary conditions be given, for use in CFD computations. Key 
measurements needed are three-component velocity profiles (mean and fluctuating), temperature profiles, surface 
pressure coefficient, and surface skin friction coefficient. 

The following final summary represents a “balanced plan” for turbulence modeling. These items are 
discussed more fully below: 

1 . Mine old experiments 

• Validate old data 

• Quantify uncertainty 

• Correlate data to particular physical phenomena 

2. Develop advanced instrumentation 

• For both wind tunnel and flight? (no consensus here) 

3. Develop effort to assess/screen existing models 

4. Assess, improve, and document numerical implementation 

5. Collaboration needed 

• Funding commitment required 
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6. Continue model development targeted to relevant flows 

• Maintain relevance to whatever collaborative plan evolves 

After assessing old datasets, new experiments may be called for to supplement, fill in, or replace where 
needed. In new experiments, an emphasis should be placed on unit problems, particularly for obtaining 3-D data. 
All old and new (proposed) experiments should be assessed as to their relevance to engineering challenges that arise 
out of this workshop. It is important to have experimental efforts both in flight as well as in wind tunnels. Data for 
validation should be organized and collected into a national or international database, similar to the existing 
European ERCOFTAC consortium effort. All experimental data should carefully provide CFD with boundary 
conditions (e.g., actual wing shape in flight, transition location, etc.). 

In the effort to assess/screen existing models, models need to be classified in some way. How should the 
assessment proceed? Should there be collaboration around multiple codes? It is probably best to have more than 
one code, but probably more than three would be unmanageable. Opening up this type of effort to too many codes 
has not worked in Europe in the past. Does the validation effort fall under NASA s domain? It might be helpful to 
form a sub-committee to devise a strategy for performing the assessment. Grid resolution and quality issues must be 
included in any study, and Navier-Stokes codes should be employed (although boundary layer codes can be useful 
both to calibrate the models and to serve as a check on the more complex implementation in the Navier-Stokes 
solvers). The validation/assessment should have a balance between old models and new ones. As a part of the 
numerical implementation assessment, consistency among different implementations of the same model should be a 
goal. Model and implementation robustness should be assessed. 

In any collaborative effort, the “right” people need to be involved. The current participants may not 
adequately represent certain segments of the aerospace industry. Other areas might have different priorities than 
those determined in the current venue. For example, the collaboration of more people working in the area of 
propulsion may be needed. As models continue to be refined and developed, any collaborative effort needs to make 
sure not to suppress new ideas that come from people outside of the “group.” Also, new models may arise both 
from experimental data as well as from mathematics and theory. 

3. Group II — Summary Findings and Recommendations. Facilitators: J.L. Thomas (LaRC) and R.A. 
Wahls (LaRC) 

3.1. Importance of the Modeling of Turbulence. The general importance of the modeling of turbulence 
is indicated below for transport aircraft vehicles and then for more general vehicle types. The key engineering 
prediction needs in the current CFD environment for transport aircraft are: 

• Reynolds Number Effects on Separation 

• Control Surface Effectiveness 

• 3-D High-Lift 

These areas are discussed more fully below. 
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3.1.1. Reynolds number effects on separation. This area is important for cruise performance, since 
increased performance is tied to designs that delay separation to higher Mach numbers or higher angles of attack. 
Transonic wings are typically designed with a mid-chord shock position and an aft-loaded section at cruise 
conditions. As the angle of attack is increased beyond cruise, the shock position moves aft, with a corresponding 
nonlinear increase in lift; a separation bubble with reattachment downstream occurs at shock Mach numbers on the 
order of 1.3. Correspondingly as the angle of attack is increased, the trailing edge separation moves forward, 
leading to a decrease in lift coefficient. These two effects are compensatory and quite sensitive to Reynolds number 
at transonic speeds. Many times, a definite change in lift curve slope occurs at the onset of separation; the 
maximum lift (and a positive pitching moment increment) occurs as the separation from the trailing edge reaches the 
separation at the shock. This sensitivity to Reynolds number is believed by some to be the root cause of the 
discrepancies between the lift levels at buffet onset in flight with the lift levels generated either by high Reynolds 
number testing or computation (see presentations by A. Kumar of NASA and M. Anderson of Boeing at this 
workshop, Appendix C). For example, with the MD-11 airframe, both wind tunnel and CFD show a much more 
pronounced break in the lift curve slope at the onset of trailing edge separation, with a consequent loss of lift in 
comparison to the flight test results at angles of attack near the observed flight buffet onset. F. Lynch of Boeing- 
retired observed that wing-body calculations on the MD- 1 1 with the Johnson-King model agreed with the trend from 
flight and differed from the SA and Menter’s k-omega shear-stress transport (SST) turbulence model results; the 
Johnson-King model has not been extended to full configurations and, thus, results for the fall configuration are not 
available. The MD-11 experience is not universal, however, since some comparisons of separated flow with both 
flight test and ground-based experiments are quite good. There is general agreement that modeling of 3-D separated 
flows is a major area of uncertainty. Reynolds number effects on separation are also especially important to 
airframe-propulsion integration problems, typified by juncture/comer regions (wing-nacelle-pylon intersections) in 
adverse pressure gradient and shock-induced separations. In general, there is an adverse Reynolds number effect 
(i.e., a decrease in effectiveness with a Reynolds number increase) associated with these flows; this effect limits the 
design tradeoffs than can be made at lower Reynolds numbers (~6 million) in transonic ground-based facilities. 

3.1.2. Control surface effectiveness. This is an area in which heavy reliance is made on ground-based 
testing because of the lack of confidence in CFD. The biggest deficiency is in 3-D applications with significant 
spanwise flow for which adverse Reynolds number effects occur. For example, F. Lynch of Boeing-retired cited an 
adverse Reynolds number effect as regards outboard aileron effectiveness on the DC- 10 wing, whereas favorable 
effects were cited for the MD-11 wing; tuft observations indicated a significantly greater spanwise component of 
flow on the DC- 10 wing. The adverse Reynolds number effects have not been encountered in 2-D. 

3.1.3. 3-D high-lift. This area is the ultimate prediction challenge because of the many flow physics 
issues involved that are sensitive to Reynolds number effects, including trailing and leading edge separations, 
confluent boundary layers and wakes, off-body separation, comer/juncture flows, and strake/chine flows. This is an 
area where extensive reliance on experiments is necessary, largely because of the lack of confidence in predicting 
separated flows. Separation on the flap is the most important driver at approach conditions, especially for advanced 
three-component high-lift configurations. Chines/strakes on the nacelle create vortices, which interact beneficially 
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with the upper surface viscous flow to control separation at high-lift conditions; these devices are generally 
determined through cut-and-try parameter variations in ground-based and flight tests. These chines and the wing- 
nacelle-pylon integration are the most important drivers to determining the maximum lift. 

The key engineering prediction needs for multiple vehicle types (i.e., military, rotorcraft, reusable launch 
vehicles, etc.) are summarized as follows: 

• Vortex Flow Breakdown and Interactions 

• Buffet 

• Active Flow Control 

• Store Separation 

• Manuever-induced Unsteadiness (Time Lags and Hysteresis) 

• Jet Impingement and Ground Interactions 

• Ducts (including Unsteady Separation) 

• Cavities 

• Rotor Blades (Turbomachinery, Helicopters) 

• High Lift 

• Transitional Flows 

• High Freestream Turbulence 

• Wake Interactions 

• Shock Boundary Layer Interactions 

• Wall Heating (Heat Transfer) 

We do not discuss these areas in detail. The above engineering prediction needs can be translated into a 

general set of flow physics issues, which cut across vehicle lines. A partial list is below. These issues could form 
the basis for a framework to classify existing experiments or advocate for new key experiments to be conducted. 

• Separation onset, progression, and reattachment, including a range of onset conditions, pressure 
gradients, crossflows, Reynolds number, and shock strength variations. The three types of separation 
typically encountered are geometry-driven (backward-facing steps), adverse pressure gradient (smooth 
surface), and shock-induced separations. Topologies of open and closed separation should be 
considered. 

• Transient evolution of transonic separated flows, including control surface deflections, Reynolds 
number, shock variations, and comer/juncture flows 

• Vortical flows, especially vortex breakdown, including stability drivers for different modes of 
breakdown, and the impact of unit Reynolds numbers 
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• Passive and active flow control devices, such as vortex generators and zero mass flow (synthetic) jets, 
including detailed data for turbulence model enhancements and the development of global, rather than 
local, models 

• Mixing layers, including merging boundary layers and wakes from the main element, flap, and slat, 
with adverse pressure gradient and Reynolds number effects 

• Curvature effects, especially recovery from curvature 

• Transition prediction and control, including trip and roughness calculations for correlation of wind 
tunnel to flight and for lower Reynolds vehicles, such as uninhabited air vehicles (UAV) 

3.2. Assessment of Current Methods. From the standpoint of vehicle prediction needs, it is clear the 
calculations involve a multitude of fluid interactions and it is often the weakest link in the elements of the overall 
process (i.e., geometry modeling, numerical method, turbulence modeling) that determines the success of the 
calculation. R. Cosner of Boeing cited the F- 1 8 wing drop phenomena as an example of the interaction of the 
various elements. The wing drop problem is a flight control difficulty that was only uncovered during flight tests; 
such surprises have a significant negative impact on program schedules and cost. The problem was addressed with 
CFD, ground-based experiments, and flight tests. The anomaly that caused the problem — an abrupt change in the 
lift curve slope at transonic speeds before maximum lift — was noted in wind tunnel tests before flight. However, it 
was judged to be a problem that would disappear at flight Reynolds numbers. The phenomena turned out to be 
Reynolds number insensitive. Initial CFD computations showed only a fraction of the lift loss that the wind tunnels 
showed and the turbulence model was the chief suspect. However, doubling the grid produced results that agreed 
much more satisfactorily with ground-based tests. At the end of the study, the program managers indicated that 
CFD gave results to within wind tunnel accuracies through the entire wing drop phenomena. This particular 
experience is not uncommon; turbulence model inadequacies are often blamed in practice for insufficient resolution 
of three-dimensional computations. 

In general, as noted in the talks presented at this workshop and elsewhere, the current capability of second 
moment closures in applications are no better than the simpler, less expensive models such as SA or SST for many 
of the key engineering predictions. For instance, extensive two-dimensional high-lift development and validation 
studies conducted through cooperative efforts by Boeing and NASA have shown that the EASM and either SA or 
SST turbulence models give quite similar results. Neither predicts maximum lift very closely, although this may be 
due to the influence of sidewall boundary layers. Studies to simulate the sidewall boundary layer should be 
completed in order to make a more definitive assessment of the turbulence model capabilities for this particular 
flow. However, there is growing confidence in the ability of Reynolds stress models to be applicable to general 
situations for which no experimental information exists. As an example, systematic comparisons of EASM 
additions to baseline models (k-epsilon or k-omega) have been recently conducted by A. Johansson of KTH [4] for a 
series of shock-induced turbulent separations for M = 5. The results show a quantum increase of accuracy over the 
baseline methods with the EASM additions, in both the length of separation and the variation of separation with 
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shock strength. In these computations, much of the improved result is attributed to improvements in the near wall 
asymptotic behavior gained through EASM. The comment was made that these flows could also probably be 
accurately computed with the SST or SA models, in which case the improvement is attributed to the variable eddy 
viscosity coefficient term rather than the nonlinear terms. Additional examples were cited of engineers at European 
car companies routinely using EASM-type models in calculations using tens of millions of grid points with notable 
improvements over linear eddy viscosity models (LEVM) for separated flows. 

As encouraging as the results using these second-moment closure methods have been, the possibility was 
discussed to circumvent the EASM class of models in favor of going directly to the RSM class. However, the 
numerical difficulties were considered so great with this class that the EASM approach should not be bypassed, 
since EASM allows many additional effects to be included rather easily into current numerical formulations. 

Until recently with the advent of direct simulations, the only way to assess turbulence models was through 
systematic comparisons with experiments. These experiments fall into two categories. The first are application 
tests, involving measurements of specific aerospace configurations, such as wing-bodies or multi-element wings, 
which are useful to the practicing engineer as a basis for the acceptance/verification of methods. These tests are 
generally not appropriate for turbulence model development, since the measurements are usually limited to surface 
measurements, such as pressures, skin friction, and oil flows. The second type are unit problem experiments, 
intended to be representative of a limited number of specific types of flow physics issues encountered in application, 
that involved more detailed measurements such as velocity profiles or turbulent shear stresses. These experiments 
can supply global or local information useful to the development of turbulent models, as in integral method 
development or simpler half- or one-equation models. For example, the Bradshaw structural coefficient of the ratio 
of shear stress to kinetic energy has been observed to hold in many flows and is used by the lag entrainment integral 
method and the Johnson-King model. These unit problem experiments can be quite expensive to conduct, but serve 
as definitive benchmark tests of turbulence models. 

For the second-moment closure-type models, experiments are not used for development, but only for 
verification/assessment, since most of the modeling is done for homogeneous flows or low Reynolds number flows. 
In this respect, the last ten years have seen an increased usage of information from direct simulations in the 
development of these methods; this trend should increase proportionally to the computational capability available 
for direct simulations. 

The turbulence modeling issues are listed below followed by a discussion of the current capabilities and 
limitations in these areas. 

• Separation and Post-separation 

• 3-D Effects 

• Unsteadiness 

• Length-scale Equation 

• Role of Curvature 
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• Reynolds Number Scaling 

• Vortical Flows 

3.2.1. Separation and post-separation. Many of the current models seem to do reasonably well for the 
few two-dimensional test cases available; this has been achieved in the models through a variable eddy viscosity 
coefficient in the formulation, which has the effect of reducing shear stress levels at separation, which tends to 
improve the correlation with experiment of the models. Current models do not uniformly predict the region 
downstream of separation, including reattachment (see [5]). Although the Bachalo and Johnson axisymmetric bump 
flow [6] is generally predicted quite well in terms of pressures (and skin friction), the shear stress levels downstream 
of separation are underpredicted. Likewise, for the backward-facing step computations, the overall extent of 
separation seems to be predicted reasonably well by some models, but all models underpredict the shear stress 
levels. 

3.2.2. 3-D effects. There is insufficient experimental data for wings with strong crossflow and separation 
effects to make a definite assessment between various models. For example, for transonic high aspect ratio wings, 
no boundary measurements are available near the trailing edge, even for attached flows. Since there is thought not 
to be a universal 3-D law of the wall behavior, methods based upon such an assumption would be less accurate. 
Thus, models that depended upon this law would have a major limitation in flows with significant crossflow effects. 
However, this limitation, if it exists, is not confirmed by the comparison of calculations with experiments to date, 
including those of 3-D calculations that use wall functions either in combination with one-equation, two-equation, or 
second-moment closures. 

3.2.3. Unsteadiness. Some transonic separated flows over airfoils with boundary layer to chord ratios of 
nearly a half were cited as quite steady with no large-scale motions evident in flow visualizations. The general 
consensus is that the current models do quite well when there is a large distinction in the time scales associated with 
the turbulence and that associated with the reduced frequencies of the aircraft motion. Several recent calculations of 
the biconvex airfoil of McDevitt at transonic speeds, characterized by shock/boundary-layer interactions which 
induce alternating upper and lower surface separations, have shown results [7, 8], which show improvements using 
SA and EASM over algebraic turbulence models. As the separation zones becomes larger, as in spoiler or bluff 
body flows, the importance of incorporating DES models increases, although there is not consensus on the extent to 
which unsteady Reynolds Averaged Navier-Stokes (RANS) by itself can be pushed. F. Lynch of Boeing-retired 
noted that unsteady flow tended to occur at transonic conditions for airfoils with small upper surface curvatures, but 
that it was only a small effect for curved sections, such as encountered on modem aft-loaded transport sections. 

3.2.4. Length-scale equation. This equation is viewed as quite ad hoc, even in attached flows and the 
contribution of the current modeling deficiencies for separated flows is currently unknown. 

3.2.5. Role of curvature. As noted in the talks presented at this workshop, additional curvature terms are 
required for LEVM models, but appear to be less important for EASM. RSM can handle the effects of curvature 
without any additional modeling. The inclusion of curvature terms in 2-D multi-element high-lift computations has 
not appreciably changed the character of the results [9]. 
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3.2.6. Reynolds number scaling. From the standpoint of turbulence modeling, it is more important to 
conduct an experiment with fully established turbulence than to conduct an experiment at high Reynolds numbers; 
the computation that reproduces the fully turbulent interaction, albeit at less than flight Reynolds number, can be 
made at higher Reynolds number with confidence. Most of the Reynolds number scaling problem is associated with 
wind tunnel to flight scaling, in which transition effects are paramount. In production testing, the objectives are 
usually quite different from establishing fully turbulent flow; for instance, it is common to locate the trips to match 
the boundary layer displacement thicknesses at the shock between wind tunnel and flight. 

3.2.7. Vortical flows. A general consensus is that all of the LEVM models fail when streamwise vorticity 
is present, such as in the computation of vortex rollup and breakdown. Even though vortices are present for all 
airframe configurations, the resolution of these vortices is generally not a significant driver for airframe performance 
prediction except if they interact strongly with the flow field, such as in vortex breakdown for low aspect ratio wings 
or chine/strake vortex interactions with high-lift configurations. 

3.3. Directions for Improvement of Turbulence Models. The goals for improvement of predictions in 
separated flows are listed below: 

• Increased generality of the formulation 

• Naturally good behavior near walls 

• Curvature corrections, especially in EASM 

• Improved modeling of active and passive flow control devices 

• More extensive evaluation of Detached Eddy Simulation (DES) methods 

• Increased role of DNS/LES for prediction modeling 

To obtain funding support for the improvement along these fronts, it is necessary to advocate on the basis 
of the improved capability that is tied to the advances in turbulence modeling. For example, efficient wing designs 
can be pushed to higher cruise Mach numbers through advances in separation prediction. Simpler and cheaper and 
more effective high-lift performance can be attained through improved computations of active flow control devices. 
Much of the modeling improvement should be tied to the validation experiments needed to assess ongoing LES 
efforts and to the assessment of capabilities across the spectrum from RANS - Unsteady RANS (URANS) - DES - 
LES - DNS. There is a need to get a national consensus on a few canonical benchmark experiments. These 
experiments should be conducted jointly with modelers and carefully designed to provide irrefutable data for the 
assessment of our current capabilities. The required experiments are expensive and should capitalize on the sizeable 
investment of other government agencies in turbulence prediction and turbulence modeling. Directions for 
improved turbulence models and a discussion of additional needed experiments are presented below. 

3.3.1. Improved turbulence models. Directions for improvement of current turbulence models are in the 
three areas below: 

• Pressure-strain Modeling 

• Near-wall Modeling 
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• Length-scale Equation Modeling 

The pressure- strain modeling pertains to models based on Reynolds stress equations. The near-wall 
modeling and length-scale equation modeling difficulties, especially the latter, pertain to most models. These areas 
are discussed below. 


3.3.1.1. Pressure-strain modeling. This modeling aspect has received extensive attention because it is the 
principal modeled term in homogeneous flows and can be studied in detail; thus, it can be judged as the most mature 
of the three areas, although the modeling in the near-wall region is still in question. 

3.3. 1.2. Near- wall modeling. The near- wall modeling problem is defined as the difficulty in the 
integration of the equations to the wall. This difficulty, usually arising in computations as an observed robustness 
problem, is the most significant limitation to usage of second-moment closure models in practice. The general 
approach is to model the equation based upon a known behavior, such as the log-law behavior for attached flows. 
At separation, this law breaks down. This is an area in which DNS is expected to be used with a reasonably high 
degree of confidence for model improvement and validation. DNS simulations have been completed with a smooth 
separation from a solid wall by prescribing the outer normal velocity in the simulation. Extension of these 
simulations to higher Reynolds number would be expected to provide an excellent source of information as to 
appropriate scaling through separation. An argument in favor of the Reynolds stress methods is that the equations 
could be formulated to be entirely independent of the wall, which one could never do with one- and two-equation 
models. A. Johansson of KTH indicated recent progress has been made in deriving and testing a second-moment 
closure method for simple flows based on realizability considerations with no wall damping terms [ 10]. This 
approach is preliminary but will be tested for more complex flows in the future. The comment was made that this 
approach should be in pressure-strain models only; one still needs damping in other parts of the equations. 

3.3.1.3. Length-scale equation modeling. The length-scale equation modeling refers to the dissipation 
equation, which is a key ingredient to all of the two-equation or higher models. It is the area of greatest uncertainty. 
The modeling is highly questionable for separated flows or for any flows with disparate length scales that interact, as 
for example, a separated airfoil trailing edge region, with momentum transfer across the wake formed by the 
merging upper and lower surface fluid. A simple mixing layer is an example where the shear stress is predicted 
accurately with k-epsilon, but turbulent kinetic energy (k) is not. Thus, the eddy viscosity coefficient is consistent 
with an explicit algebraic stress model but one has to change epsilon to get k correct, clearly indicating a problem in 
the length scale equation. The length scale model changes as the models are changed from LEVM to EASM to RSM 
(i.e., production terms and diffusion terms are known slightly better with RSM since some of the terms are 
computed directly). Flows in which turbulent transport is important clearly require RSM approaches; however, the 
length-scale equation deficiency may be masking effects that should be otherwise accurate and computable. 
Advancements could be made by using DNS or spectral theory to suggest models accompanied by LES 
computations or experiments for model validation. The main concern with this area, as with all methods based on 
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DNS or LES, is the scaling of the results to flight Reynolds number associated with flight vehicles, since these 
simulations will not be practical for many years. 

3.3.2. Needed experiments. The need of measurements in unit problems was discussed at length. These 
unit problems are used to demonstrate model improvements in treating deficiencies of current models, such as in 
post-separation regions. There are only a few test cases available to assess turbulence models for flow separation 
that are not plagued by three-dimensional interference effects. These include the Bachalo and Johnson 
axisymmetric bump [6], the Driver flow [11], and the separated flow of Simpson [12]. The latter two flows are very 
difficult to simulate with Navier-Stokes computations because of boundary conditions, but could be computed using 
inverse techniques to determine an effective wall shape, such as that used by Rumsey and Gatski [13]. Two 
dimensional axisymmetric flows (delta/R « 1) are recommended as unit problems since there is less influence of 3- 
D effects, as for instance from sidewall boundary layers and secondary flows, and the data collection is easier than a 
3-D flow. These experiments could build upon the Bachalo and Johnson axisymmetric bump experiment; this 
experiment is characterized as a trailing-edge separation experiment in which the separation moves upstream to the 
shock and has been widely used to compare different models. It could be expanded using current CFD design 
methods to include shock-induced boundary layer separation or modified to induce controlled three-dimensional 
effects (through addition of sweep, suction, or a vortex generator, for example). The experiment would be simple 
enough in complexity to be computed by LES and/or DNS and, thus, serve as a validation case for such approaches 
at lower Reynolds number. 

There were no specific recommendations proposed for 3-D experiments to serve as definitive benchmark 
datasets. Difficulties arise in measuring boundary layer information at transonic speeds at reasonable Reynolds 
numbers. In this regard, no velocity or shear stress data for transonic flows with strong 3-D effects at high Reynolds 
numbers is available. The NASA trapezoidal wing, high-wing, test case is viewed as a step along the way for 3-D 
high-lift validation, but lacks the key element of nacelle-pylon-chine integrations encountered on realistic airframes. 
Also, the flow is quite sensitive to the component rigging of the high-lift configuration and small adjustments can 
accentuate or mask certain intended flow interactions. 

Adaptive turbulence models (also called zonal turbulence models)— using various turbulence models as 
they are appropriate to the local physics— are a possible way to circumvent the lack of a general purpose turbulence 
model applicable to all flows. Such techniques are not widely used largely because there is no a priori knowledge of 
the capabilities of a given turbulence model in a given situation. Establishment of a basis to determine the capability 
of a given model in a calculation would be quite useful. Methods to determine the error associated with the 
discretization of a given set of partial differential equations and use that as a basis for adapting the grids to attain a 
specified error tolerance are now being pursued. Also, guidelines/standards to ensure sufficiently accurate 3-D 
computations are now emerging, as noted in the recent ERCOFTAC [3] referenced above. However, capabilities for 
determining the accuracy of complex-geometry separated flow computations are just in their infancy. Moreover, for 
complex flow applications at high Reynolds number, there seems to be nothing on the horizon to account for the 
physical error embodied in the solution of a given set of turbulence models. 
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3.4. Conclusions and Recommendations. The general recommendations for this group focused on 
directions for improving turbulence models. The highest priority should be given to attacking the length-scale 
equation, since it is a key ingredient to all of the two-equation or higher models. Advancements could be made by 
using direct simulations and/or spectral theory to suggest models and then testing these models by comparing to 
DNS/LES computations or model validation experiments. The second priority should be in developing improved 
near- wall behavior, since the stiffness issues associated with this region pose a significant limitation to usage of 
second-moment closure models in practice. In this area, direct simulations could be used with a high degree of 
confidence for both model suggestion and validation. A principal concern with this approach is the ability to scale 
the results to flight Reynolds number. 

4. Group III - Summary Findings and Recommendations. Facilitators: W.L. Sellers (LaRC) and R. 
Rubinstein (LaRC) 

4.1. Importance of the Modeling of Turbulence 

4.1.1. Vortical flows. Vortical flows are ubiquitous in aerodynamics. Flows around vortex generators are 
common examples. An interesting non-aerodynamic application of current interest is the hydrocycle centrifuge for 
removing impurities. Unsteady vortex bursting on fighter aircraft was identified as a source of fatigue damage 
through flow-structural coupling. 

4.1.2. Separated/time-accurate flows. High Reynolds number separated flows occur in buffet onset. 
The claim that all models fail to predict the dependence of the lift coefficient on angle of attack helped motivate the 
present workshop. These flows also occur in high-lift airfoil configurations and in the evaluation of control surface 
effectiveness. 

4.1.3. Juncture and mixing flows. Aerodynamic juncture flows include the flow around the wing-body 
nacelle strut, and the flows around pods and blisters, which are particularly important for military aircraft. 

Mixing flows occur in confluent boundary layers, wake interactions with solid bodies in high-lift 
configurations, and in turbulent wake impingement in turbomachinery. Potential difficulties exist because the correct 
modeling of mixing of disparate turbulent flows remains somewhat obscure. If we think of a homogeneous region of 
turbulent flow as generated by forcing at some integral scale, then when two turbulent flows mix, it is not clear how 
the integral scale of the result should be determined. In unit problems like the penetration of turbulence generated by 
an oscillating grid into quiescent fluid, current models are frequently found inadequate. 

4.2. Assessment of Current Methods 

4.2.1. Vortical flows. The first general observation made was that vortex flows could pose significant 
numerical issues that are independent of any turbulent modeling consideration. Insufficient grid resolution alone will 
diffuse small tight vortices. There is considerable numerical research on adaptive gridding methods to identify and 
resolve vortices. The implementation of such methods on unstructured grids is especially important, and was 
addressed in some papers at the 200 1 AI AA January meeting. 
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Otherwise, as noted earlier, vorticity-dominated flows present multitudinous problems for turbulence 
models. Eddy viscosity models (including the quadratically nonlinear models) incorrectly predict a rigid rotation 
profile of swirl velocity inside vortices. Cubic nonlinearity was suggested as a solution to this problem. 

The Reynolds stress model (RSM) reduces the turbulent kinetic energy inside the vortex, correctly 
indicating the suppression of turbulence in the vortex core. This effect is not predicted by lower order models, even 
by lower order models with swirl corrections. Examples were given by Dr. Kim of FLUENT during the 2001 AIAA 
January meeting. 

The flow around a 6: 1 prolate spheroid provides an example of a vortical flow, since vortex interactions 
can occur, leading to unsymmetric forces on the body. Unsteady separation of these flows is important. There is a 
European effort in progress to calculate this flow. Preliminary findings are that the eddy viscosity models (EVM) are 
generally inadequate, but that the RSM is at least qualitatively correct. 

4.2.2. Separated/time-accurate flows. Turbulent separated flows have proven to be difficult to 
compute. Earlier discussion of mixing and juncture flows, which are otherwise amenable to all current turbulence 
models, showed that separated cases of even these flows are hard to compute accurately. 

There was general agreement that no model, even RSM, can predict recovery of turbulence downstream of 
separation. Associated with this fact are the observations that typically the models overpredict the turbulent time- 
scale, and that turbulence production takes too long to re-establish itself. It was noted more generally, that in many 
problems involving the recovery of turbulence, all models are inadequate. Recovery following shock interactions 
also belongs to this class. Very large-scale structures and long time scales characterize these flows. Consequently, 
even in DNS, it is necessary to integrate for a very long time to resolve the post-separation region. 

Another motivation for this workshop had been the suggestion that there is a Reynolds number effect that is 
missing from all current models. Thus, if a code is validated against ground data, it may fail for flight conditions just 
because the Reynolds number is larger. After very careful discussion of this point, it was agreed that Reynolds 
number does certainly influence quantities like the thickness of boundary layers and the location of transition. 
However, existing models are properly sensitive to such effects. It was agreed that no other Reynolds number effect 
exists. 

A question was raised about the accuracy required in calculations of lift and drag coefficients. Accuracy to 
the order of 1-2% seemed to be desired. However, none of the nominally universal constants of turbulence theory, 
such as the Kolmogorov constant, or the eddy-viscosity constant C\x is known to within less than 10%. It is possible 
that the calculations are not very sensitive to the modeling constants, or that the model constants are empirically 
adjusted. Nevertheless, this question deserves serious consideration. 

4.2.3. Juncture and mixing flows 

4.2.3.I. Juncture flows. The main difficulty in predicting juncture flows is attributed to flow separation, 
particularly in 3-D flows. In the wing-body nacelle strut, local separation occurs, leading to a 1-2% drag increase. 
Two-equation models predict incorrect primary separation. The SST model gives better separation predictions, but 
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cannot predict secondary flow; however, secondary flow is not always important. For these flows, the RSM is 
adequate, and the EASM should work. 

In 2-D separated flows of this class, eddy viscosity models such as the SA and SST are adequate, and in 
fact, no models are better. But in these cases, normal stresses are not important: prediction of normal stresses always 
requires higher-order models like the EASM or RSM. 

An entirely different application of this type of flow is in icing calculations. It was noted that heat transfer 
is often underpredicted in these flows, leading to underprediction of ice growth. 

Juncture flows often exhibit unsteady, chaotic, large-scale vortices. These vortices contribute to mixing and 
cannot be ignored. Nevertheless, there was speculation than any Reynolds-averaged model may suppress this type of 
effect. Phase-averaged models are a possible alternative. 

4.2.3. 2. Mixing flows. In the computation of confluent boundary layers, all current models (SA, SST, k-co, 
and EASM) provide satisfactory results for both the mean flow and turbulence quantities as long as the flow remains 
attached. The difficulties encountered with separated mixing flows must be attributed to separation rather than to 
mixing. 

Similarly, for wake flows, including the wake from the main element of multi-element airfoils, all models 
provide satisfactory predictions in the near-wake region. In high-lift computations more generally, it is found that 
for low angle of attack, the unsteady slat flow is predicted well. In wind tunnels, this flow is transitional, and the 
transition prediction poses difficulties for models; however, in flight conditions, transition is not a problem, and 
models agree well with flight data. 

4.3. Directions for Improvement of Turbulence Models 

4.3.1. Vortical flows. Swirl and vorticity corrections to two-equation models were discussed 
extensively, but ultimately rejected as less satisfactory than EASM and RSM. 

4.3.2. Separated/time-accurate flows. Further investigation of the failure of models to predict the re- 
establishment of turbulence is needed. Unfortunately, it was also agreed that no one really knew what specific 
weakness of current models is responsible for the problem. 

4.3.3. Juncture and mixing flows. The difficulties posed by unsteady vortices raised the general issue 
that turbulence models are always more successful in predicting strain-dominated flows than vorticity-dominated 
flows. Some vortex flows are actually predicted better by purely inviscid methods. We need some way to identify 
these problems in advance, perhaps through a discriminant function. 

4.4. Conclusions and Recommendations 

4.4.1. Vortical flows. In general, RSM and EASM are satisfactory solutions to this class of problems. 
However, the numerical issues noted above remain. Continued progress on the numerical side will be important. 
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4.4.2. Separated/time-accurate flows. The existence of a “Reynolds number scaling” problem, 
understood as an intrinsic failure of turbulence models due only to high Reynolds number, was emphatically and 
unanimously rejected. 

4.4.3. Juncture and mixing flows. In the absence of transition and separation, current models are 
adequate for aerodynamic mixing flows. 

A possible limitation of Reynolds-averaged models was discussed, although there was not general 
agreement about this possibility. 

Better models for vorticity-dominated flows, and methods for distinguishing strain-dominated from 
vorticity-dominated flows are needed. 

A list of problems in which large discrepancies exist between CFD and data should be compiled. 

5. Conclusions. This section is divided into two parts: an assessment of current methods, and 

recommendations for future development. These represent a collation of the key points and final summaries from 
the three groups. More detail regarding these assessments and recommendations can be found in the groups’ 
respective sections above. 

5.1. Assessment of Current Methods. 

1. Key areas, or flow categories, in the prediction of aerospace vehicle performance for which turbulence 
modeling improvements are particularly critical: 

a. Separation, particularly separation onset and progression in three-dimensional flows 

b. Vortex flows, particularly vortex breakdown and component interactions 

c. Jets and mixing layers 

d. Unsteady flows 

2. No current model can predict post-separation (recovery) physics, and most models cannot predict post- 
curvature (recovery) physics. 

3. There is no evidence that Reynolds number scaling is an issue in turbulence models in the sense that if a 
model does well at low Reynolds number, it should have no trouble at high Reynolds number. However, 
ground-to-flight scaling remains an important issue because of transition and other effects. For a given 
configuration, the severity of adverse pressure gradients and the degree of spanwise flow increases with 
increasing Reynolds number. Thus, it may not be sufficient to validate turbulence models only at lower 
Reynolds numbers, and high Reynolds number experiments remain a critical part of the validation process. 

4. Turbulence models are often blamed for problems that are caused by insufficient grid resolution, 
unaccounted-for transition effects, or geometry issues (such as aeroelastics). Isolation of turbulence model 
deficiencies requires a critical account of all of these effects. 

5. The current capability of second moment closures in aerospace applications are no better than the best of the 
simpler, less expensive models such as SA or SST for a large number of key engineering predictions typically 
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dominated by thin shear layers. Hdwever, there is growing confidence in the ability of Reynolds stress models 
to compute not only this class of flows but also more complex turbulent flow classes. 

6. There may be fundamental limitations in the ability of RANS approaches to predict critical industry needs for 
separation onset and progression at flight conditions. However, given that computer limitations force the use 
of RANS methods for realistic configurations at high Reynolds number in the foreseeable future, the ASM 
and RSM models with much additional development currently offer the best hope. 

5.2. Recommendations for Future Developments. 

1. Establish a national/intemational database and standards for turbulence modeling assessment 

2. Focus modeling efforts on length-scale and near-wall modeling, with emphasis on second-moment closure 
(and derivative) models 

The workshop participants strongly advocated either establishing a new database of high quality, well- 
documented experimental data or else tying into the European database (ERCOFTAC) already underway. The 
database should include “unit problems" that are simple and isolate specific aspects of turbulence. A thorough and 
exhaustive search and categorization of existing experimental databases should be conducted prior to creating new 
ones. New experiments, for instance three-dimensional flows with separation and crossflow effects, especially at 
high Reynolds numbers, can then be performed to either upgrade or fill in missing datasets. Continued development 
of advanced instrumentation would benefit future experiments by improving the quality of the data, or allowing 
quantities to be measured that currently cannot be easily obtained. 

There are still many uncertainties regarding the precise failings and capabilities of existing turbulence 
models due to non-uniformity in coding and evaluation practices. The development of a standard process for 
assessment would be beneficial (see, e.g., [14, 15]). The workshop participants also identified the fact that numerical 
methods for model implementation and model-dependent mesh resolution requirements are relevant issues in the 
quest to improve turbulence models. Therefore, a part of the standard should include assessment and documentation 
of the mesh requirements and numerical implementation. The latter issue is important because of current non- 
umformity in programming practices, which sometimes results in ostensibly the same model performing differently 
in different codes. 

Continued turbulence modeling research and development is important, particularly in the areas of second- 
moment closure models such as RSM and EASM. Future development efforts should maintain relevance to the 
specific needs of the aerospace community. A significantly increased role for DNS and LES/DES for prediction 
modeling is possible because of emerging simulations with separated flows. Although still limited to low Reynolds 
numbers, DNS and LES should be able to serve as a foundation for the development of improved RANS models for 
the length-scale equation and the near wall treatment, because they can provide information on quantities which 
have a crucial role in turbulence modeling, but which are difficult to measure in laboratory experiments. The 
length-scale equation is a major deficiency and area of uncertainly with most current models from two-equation to 
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second-moment closure models. The modeling improvements could be used in the practical computations of flight 
vehicles at Reynolds number beyond the range of simulations. 

Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent viable 
approaches for analysis on a limited basis, and should remain a part of any balanced effort in turbulence modeling. 
These methods will be required for accurate predictions of highly unsteady turbulent flows in which scales of the 
Reynolds-averaged velocities and turbulent velocities are not widely separated. 
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over the complete flight envelope, thus allowing greater 
of computational tools in the aircraft design process 
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3/25/98 - Cosner - Boeing PW CFD Cost and Cycle Time Improvement 
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NASA Flight Prediction Roadmap (Near-Term) 
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Turbulence Modeling Workshop 

Expectations 
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figures or Merits to evaluate turbulence models, e.g. a unit 
problem, complex problem, robustness of the model, integration 
with computational codes, ease of use by a designer, etc. 



Workshop Structure and Outcome 
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Assessment of today's status from Subcommittee perspective. 



NASA ASTAC Airframe Systems Subcommittee Perspective 



50 


Sub-Com. ^ Division Director 
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Aircraft capacity balance (regional / business aircraft vs. commercial transports). 



NASA ASTAC Airframe Systems Subcommittee Perspective 
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'Current CFD (turbulence modeling) capabilities do not permit the reliable prediction 
of separation onset / progression characteristics. Efforts currently underway to 
address this deficiency, which is a critical path requirement if desired design cycle 
time reductions are to be achieved, need to be accelerated. Major advancements are 

needed.' (Letter to ASTAC Chair, 21 March 2000) 
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NASA ASTAC Airframe Systems Subcommittee Perspective 
12 January 2001 ICASE / NASA Turbulence Modeling Workshop 
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Assess prospects for Reynolds Averaged Navier Stokes (RANS) turbulence 
modeling vs. alternate approaches (DNS, LES, etc.) 
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NASA / Boeing CFD Buffet Onset Team has made significant progress. Much yet 

remains to be done. 

• Outstanding Subcommittee concerns should be pursued vigorously within 
proposed 21st Century Vehicle Technology Program. 
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We must work together. 
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Outlook. 




Ultimate Expectations from CFD 


c 

0 

D 

CT 0 


o o 
— c /: 
0 to co 
c *- "o 

0 ^ ^ 

1 * T" L. 

O m 0 

■— 0/ c\ 

TJ£ r 

0 ■*- c 

L C 3 
Q_ - C= 

<f) — 

-*— » »— -£— 

o o o 

0 i= ro 

t in ^ 
0 . 


co g 

s/ O 

■■■ ■ t *v 


2 o = 


0 b 
Q. . 2 
E th 
O^D) 
o c c 

CO CO - 0 

' 0 ° 

Q — CD 
LL CD “ 

^11 

S’s.i’ 

O J_ X3 

T (D r 5 
*■* >* 13 
C CD D 
0) T A 

> ^ 2 
0 b 
CD m 
(/) D 1 - 
“ C CD 

lol 

hCOS 

I I I 


r •» ® 

■§> o 

e e? 

O 0 O 


C i 
— "D 
0 0 
0 1 I 

E ° 
.E ro 


>> t 

0 . & 
£ a> c 
0 E o 
Jr; "O 

>< 0 C 

0 _£ 0 

co ^ Q. 
— = 0 
(/) CD "D 
O K — 
0 o 

iS 3 C 
0 CO O 

o “ o 


-5 

• ■ o 
0 a 

o S" 

O 'M 

2 x 

o • 

4— O 

O -D 
_C v 
O *5. 

CO Q. 
O CD 

ro g 

0 o 


O 0 

o ^ 

c" CD 

CD ~ 
c* C 


S 5 t 

0 >* 


M 8 O 

1 si 

C (D 5 


0 0 
c Q 

8 o 


f t=I 

«_ -c 
0 o 5 

_C o Q- 
.ti r 0 

LU .2 CO 


0 

c 

. ■ o 

0 --F 

£ o 
3 > 

0 2 

0 0 

4— “ 

=> ° 

> a. 

_o w 

jco 

"O 0 
C 0 
CO ^ 
0 o 5 

0 D) 

o 5 

> * 

0 

O .E 

’E o 

CO _0 

E “0 

Q 

■a cd 
o z 


59 


RANS modeling of free vortices is very poor (diffusive is far off). 

For VG’s, treated collectively (“manipulated TBL”), it’s non-existent. 








Wing Twist in Cryogenic Tunnel 
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Small Flow Features Have Much Control 
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(who needs adaptive grids?) 





The Nature of RANS Turbulence Modeling 
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LES is decades away. 

Not to mention Quasi-DNS... (“LES” with viscous sublayer resolved) 



Airfoil Testing 
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suction 0 /Mifiwg Jiang, Long Beach, 2000-0510 









Effect of Karman Constant on Skin Friction 
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The Challenge of 2D Boundary Layers 
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Unsteady, Three-Dimensional, High Reynolds. 

Truly more challenging? 
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The same for “high Reynolds number”: no mystery there. 

— K matters! 



Simple versus Elaborate RANS Models 
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but it could be only in boundary layers {next slide). 


Detached-Eddy Simulation 

The words “demonstrable progress” in the invitation were mine. 
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stall/spin, spoilers, landing gears, cavities, noise.. 



The Resolved Solution in Different Approaches 
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Cylinder, laminar separation, 1.e5 il^lliits g 



Circular Cylinder: Grid Convergence?! 

Laminar Separation, /?e®tsl0 5 
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URANS 
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LM Aeronautics Transition 
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M Aero - Fort Wo 



-LM Aeronautics Lines of Business: 
acAir, Air Mobility, & Special Projects 
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Turbulence Modeling is a Major Source of 
Error in CFD Simulations 
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Turbulence Modeling Accuracy Requirements: 





Fighter Aircraft Turbulent Environments 
Complex Geometries - Complex Flows 
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Advanced Fighter Aircraft and UCAV 
Internal Flows - Radical Designs 
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Axial Location (X / L) 
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Turbulence Modeling Practice in 
Aerospace Industry 
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• Experience base and expectations with current models make 
change difficult. 



Shortcomings of Two Equation Models 
for Aerospace Applications 
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S ear term implementation in the industrial environment, models that 
ase computational requirements by more than a factor of 2 or 3 
J likely be useful only for special cases 



Turbulence Model Accuracy for Simple 
Flows is a Major Source of Concern 
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ne Test-to-CFD Comparison 

quation model underoredicts losses in 



Axial Location (X / L) 



for Test-to-CFD Comparison 

equation model closely matches performance with flow control. 



Axial Location (X / L) 


wb- Improved RANS Models Can Account 
< for Complex Turbulent Phenomena 

Non-equilibrium two equation 
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strophy model 
\ • Complex 
\ Not mature 
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Large Eddy Simulation Necessary for 
Some Aerospace Applications 
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Simulation of Weapons Bay for Structural jjf 
Loads Requires Unsteady Analysis 

Symmetry plane Mach contours show 

ssures needed for structural resolved unsteady structures in flowfield 
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Pod Flow Field Characterization 
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f-Pod Wake Acoustic Signature 



-OW Speed Wind Tunnel 91 118 0.292 

\ Water Tunnel PIV 0.65 0.92 0.268 

fii^ve Accurate CFD:LES 150 223 0.248 


















LES Modeling Development 


=3 o 


CO t± 

W Q_ CO 

o 2 0 

•^3 Q. 0 

_D CO >* 

O ID CO 

w co 00 

S T3 W 

□_ q Q. 

m -Q 
Jr i- 0 

CD O > 

0 CO o 

CO -Q _r 

= CD ^ 

13 CO 42 

*+- m 

.£ -= 0 

■2 ■£ o 

1 is 

O Q. CO 

: ofl) 

1/5 0 r^v 

O > O) 

O 0 .E 

• "0 0 

CO "0 
111 ° 
_i E 


% 


Page 23 


Outline 
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Experimental Accuracy Considerations 
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Turbulence levels 
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Directions for Improved Turbulence 

Modeling 
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Validation and Documentation 
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• Requires long term commitment of budget and resources to 
maintain database 
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Wind Tunnel and DNS Data 
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Model Implementation 
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Conclusion 
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TURBULENCE MODELING: 
A NASA PERSPECTIVE 


T. B. Gatski 

Computational Modeling & Simulation Branch 
NASA Langley Research Center 
Hampton, VA 23681, USA 


Turbulence Modeling Workshop 

January 12 — 13, 2000 
Reno, Nevada 
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OBJECTIVE 

• Identify deficiencies in predictive capabilities 

• Develop improved/new models 

* Capable of more accurate predictions of flows currently of interest 

* Capable of solving more complex flows which become of interest 


APPROACH 
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HIERARCHY OF SOLUTION METHODS 



• RANS approach currently most common methodology for calculating 
turbulent (aerodynamic) flow fields 

• RANS formulations inherently susceptible to the closure problem 

* As statistical moment equations are derived, higher-order moments 
appear that require closure 

• Within RANS framework, level of sophistication used in developing 
the turbulent closures varies widely 

• A separate though related issue is the numerical solution of the closed 
set of equations needed in RANS formulation 

* CFD issues associated with accuracy and efficiency of solution of 
RANS equations depend on closure level 
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TURBULENCE MODELING FOCUS 

• Single-point (space and time) correlations 

* Linear eddy viscosity model (LEVM) 

* Nonlinear eddy viscosity model (NLEVM) 

§ Algebraic stress model (E)ASM 

* Second moment closure (SMC) 

§ Reynolds stress model (RSM) 

SINGLE-POINT CLOSURES 

COST, 

COMPLEXITY 

REYNOLDS STRESS 

/ EQUILIBRIUM FOR ANISOTROPY 
/ SIMPLIFIED DIFFUSION TENSOR 

ALGEBRAIC STRESS 

SOLUTION OF TENSORIAL EQUATION 

u 

FORMAL 

EXPLICIT ALGEBRAIC STRESS 

“PHYSICS ', -i DYNAMIC RANGE" 

• Closure model development focused in incompressible regime 

• Mean variable density extensions used mainly in compressible flows 

* Role of compressible correlations uncertain 

* Direct numerical simulation of supersonic flow (Ames, Langley) 


JT non-linear eddy viscosity 

/ -MOST- GENERAL EXPRESSION 

LINEAR EDDY-VISCOSITY 
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T' * 5‘ ' 

h — 13 — -3- 

ij 2 K 3 


* fit is the turbulent eddy viscosity 

* r is a turbulent time scale 


* b t j is the Reynolds stress anisotropy tensor 

* Term proportional to the kinetic energy K is (formally) 
absorbed into the pressure term for p 

• Coupling between mean flow and turbulence is through pt 


in 
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TWO-EQUATION MODEL 

• Many two-equation models 

* K — e model 

* K — lu model 

* Shear Stress Transport model (Ames) 

* More 

ONE-EQUATION MODEL 

• Transport equation for the turbulent eddy viscosity 

• Popular among industrial users due to its ease of implementation, rel- 
atively inexpensive cost, and good performance 

HALF-EQUATION MODEL 

• Outer eddy viscosity is modified to account for the effect of streamwise 
evolution of the flow on the turbulence 

* Ordinary differential equation is solved for the streamwise 
evolution of the maximum shear stress (Ames) 

ZERO-EQUATION (ALGEBRAIC) MODEL 

• Based on Prandtl’s (1925) mixing length theory 

• Isotropic eddy viscosity 
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LIMITATIONS AND ALTERNATIVES TO LEVMs 

• LEVM are a proven tool in turbulent flow-field predictions 

* Inherent in the formulation are several deficiencies 

§ Isotropy of the eddy viscosity 

• Consequence of Boussinesq approximation which assumes 
direct proportionality between r x j and S XJ 

• Precludes prediction of turbulent secondary motions in ducts 
§ Material-frame indifference of the models 

• Consequence of sole dependence on (frame-indifferent) S r j 

• Insensitivity of turbulence to noninertial effects 

• Need dependence on rotation rate tensor W t j 

* Such defects cannot be fixed in a rigorous manner 

• Nonlinear eddy viscosity models (NLEVMs) and their subset algebraic 
stress models (AS Ms) extend the range of applicability of LEVMs 

* Replace the Boussinesq approximation r x j = r l j(S lJx r) with 

r l:j Aj ( Sfci , fb/c/ ) ) 

* Need for transport equations for the characteristic turbulent scales 

* NLEVM and EASM formulations based on a two-equation closure 

§ Any other (lower) closure (zero-, half-, one-, or two-equation) 
could be connected to a NLEVM or ASM. 
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NONLINEAR EDDY VISCOSITY/ 
ALGEBRAIC STRESS MODELS 


7 


• Class of nonlinear eddy viscosity models that extend the one-term 
tensor representation in terms of S t j used in LEVMs to 

nj = \k6, 3 + X a' nT V 

O n=l 

* (n = 1 ... n) is a given tensor basis, n finite 

* a' n expansion coefficients which need to be determined 


• For general Reynolds stress representations, coupling to mean flow 
can be either through the direct use of T{j in momentum equation or 
through a modified form given by 


P- 


Du, 

Dt 


dp d 

dx. 


+ 


dxj 


[ft + fk) 


duj 

dx n 


* S represents nonlinear (source) terms from the 
tion 


+ 5 

tensor represent a- 


• Generally assumed in developing closures for Reynolds stresses 

* Functional dependency on the characteristic turbulent scales 

* Functional dependency on the mean velocity gradient 

bij = b lJ (Ski, Wki, r ) 

• Tensor representations basis also assumed to be functions of S and W 
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• For b(S, W, r) this basis consists of the elements 


= S 

t( 2 ) = SW - WS 
= S 2 - |{S 2 }I 
T ( 4 ) = w 2 - |{W 2 }I 
T (5) = WS 2 - S 2 W 


T (e) = w 2 s + SW 2 - §{SW 2 }I 
- wsw 2 - W 2 SW 
= sws 2 - S 2 WS 

T( 9 ) = w 2 S 2 + S 2 W 2 - |{S 2 W 2 }I 
T( 10 ) = WS 2 W 2 - W 2 S 2 W. 


• Nonlinear eddy viscosity models (NLEVMs) 

* Expansion coefficients determined from 

§ Calibrations with experimental or numerical data 
§ Some physical consistency constraints 

• Explicit algebraic stress models (EASMs) 

* Expansion coefficients derived consistent with the results of tensor 
representations from the full differential RSM 

• In both models 

* Explicit tensor representation for b is obtained 

§ Subset of full representation basis 


115 


9 


SOME EXAMPLES 

• Quadratic NLEVM (Glenn) 

* RDT (rapid distortion theory) result for rapidly rotating turbulence 
(no shear) used 

* ol 1 coefficients determined by applying realizability constraints to 
the cases of axisymmetric expansion and contraction 

* Coefficients optimized by comparison with experiment and numer- 
ical simulation 

• Quadratic EASM (Langley) 

* Extracted form from RSM with SSG pressure-strain rate correlation 
model 

• Enhancements to improve predictive capability 

* Rotational and curvature effects (CTR, Langley) 

* Wall proximity effects (CTR, Langley) 
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LIMITATIONS AND ALTERNATIVES TO NLEVM/EASMs 

• NLEVM/EASMs are now being used in many applications 

• Inherent in the formulation are two (possible) deficiencies 

* Weak equilibrium assumption (Db t j/Dt = 0) 

* Assumed form of turbulent transport and viscous diffusion model 

• Such defects can be addressed in a rigorous manner with EASMs 

* Improving frame-invariance property of formulation (Langley) 

* Modifying assumed form for turbulent transport model 

• Close linkage between the (E)ASMs and the Reynolds stress equations 

* Directly incorporate models for the pressure strain rate correlation 
and anisotropic dissipation rate 

• Unfortunately not all features of differential SMC retained 

* Relaxation effects of individual stress components precluded 

* Turbulent transport effects only partially taken into account 
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SECOND MOMENT CLOSURES (SMCs) 

• Reynolds stress model (RSM) most common SMC at this time 

* Structure based models account for dimensionality of turbulence 

• Calibrations extensively based on homogeneous flows 

* Applied to inhomogeneous flows 


• Tij = U'ili'j 


Incompr -essible Flow 


Dt H 

Dt 


Pij + Ilij + Dij 



P = 

1 IJ 


duj 

Tlk PP r jk 


dx k 


dut 

dx k 


Ilij = p' 


du[ du'A 


dxj 


+ 


Ox, 


1} 


D ij = 


d_ 

dxi 


pu'jUk +p'(u' i 5 j k + u'jSik) - (a' ik u'j + CTjkP) 


turbulent transport 




2 

= -e8, ] + 2ed, J =a’ k ^ + a' jk ^ 


viscous diffusion 

dp dP- 


2v 


dx k dx k 


* Development of improved closure models for 

§ Pressure-strain rate correlation (Langley, Glenn) 
§ Anisotropic dissipation rate (Langley) 

§ Turbulent transport term (Langley) 

§ Wall proximity effects (CTR) 


118 


Cornpressible Flow 


• pr t j = pu'-Uj (Favre variables) 


P 


Dr^ 

Dt 


pPij ~F p Hi) + pDij 


p€ij + pMij 


_ ~ _ duj _ du 

P*ij P7~ik ^ _ P7~jk 


dx k 


dx k 



pDij — 


d 


dx k 


pUiUju'l + p'^Sjk + u'jSik) - (a ’ ik u'j + a' jk u'i) 


turbulent transport 


viscous diffusion 


2 du'j du ■ 

/96fj = -p£Oij + 2 pea \j = 




pM t j = pu" 


5(7 




3x/t 


Sp'l , _—(da ik dp 

8^ 


Development of improved closure models for 


* Dilatation dissipation 

* Pressure dilatation 

* Mass flux model 

* Wall proximity effects 
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SUMMARY OF CURRENT ACTIVITIES 

Theoretical Approaches 

• Half-equation and two-equation model improvements 

* Eliminate coordinate dependence in required “switching function” 
for SST model (Ames) 

* Extend lange of applicability of turbulence models to transition 
regime (Glenn, Langley) 

* Improve predictive capability in APG and free shear flows (Ames) 

• Improve explicit algebraic stress formulations (Langley, CTR) 

* Improve frame-invariance properties of model (Langley) 

* Improve sensitivity to rotation and curvature (CTR, Langley) 

* Formulate to account for compressible terms (Langley) 

• Some selective modeling of current RSM closures 

* Turbulent transport terms (Langley) 

* Improve wall proximity effects (CTR) 

• Structure function modeling (Ames) 

• Length-scale equation from multi-point analysis (Glenn) 

• Generalized wall functions (Glenn) 
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SUMMARY OF CURRENT ACTIVITIES 

Applications 

• Variety of aerodynamic flows 

* High-lift flow field prediction (Ames, Langley) 

* Transonic buffet onset analysis (Langley) 

* Trap- wing CFD validation (Ames, Langley) 

* Transonic vortex flow over delta wing 

* Ramjet/scramjet flowpath analysis (Glenn, Langley) 

* Variety of separated flow studies 

• Variety of unit problems (Ames, Langley) 

* Wakes in pressure gradient 

* Flows with curvature 

• Deficiencies in model performance cycled back into improved model 
development 

Other RANS Modeling Areas 

• Multi-scale turbulence models 

• Non-equilibrium effects 

• Dynamic zonal modeling 

* Linkage between LES and non-stationary RANS-type closures 

• Heat transfer and reacting flow modeling 
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STANFORD UNIVERSITY 
Thermosciences Division, M.E. Dept 
Stanford, CA 94305-3030 


ICASE/LaRC Turbulence modeling workshop 

Turbulence Modeling for 
High Reynolds Number Separated Flows 

Peter Bradshaw 

bradshaw@vk.stanford.edu 

12 Jan. 2001 

(no oral presentation) 
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MY TAKE ON SEPARATED FLOWS-I 


• (I) Virtually all Reynolds-averaged turbulence models (or SGS 
models) are calibrated in flows dominated by shear layers (in suit- 
able axes) 

• (II) There is no guarantee that any model will be "universal" 
(i.e. reliable in flows very different from those used for calibration) 

• (III) (important parts of) Separated flows are not dominated by 
shear layers 

• (IV) Therefore Reynolds-averaged models cannot be guaranteed 
in separated flows. Ever. 
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MY TAKE ON SEPARATED FLOWS-II 


• The problem of predicting separated flows has two parts:- 

• (A) predicting the separation line* with a more-or-less given 
pressure distribution 

• (B) predicting the separated region with a more-or-less given 
separation line 

. . . “more or less” because there is strong upstream influence - but 
initial tests of a model could and should be divided into (A) and 

(B) 

• “line” not “point”: all real-life separated flows are three-dimensional, 
some with very strong streamwise vorticity 
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TEST CASES FOR SEPARATED FLOWS 


• Predicting reversed -flow skin friction in the 2D backstep flow is 
still a useful test. . . 

. . .because it goes as Re -0 - 5 approx. (0.5 ^ 1/2) and this should 
be reproduced by a low-Re (wall-layer) model 

• The largely-neglected flow over a cone at zero incidence is a 
severe test of models for the separated region because that is so 
far from “still air”. (Calvert, JFM 27, 273 (1967): probably only 
the base pressure is reliable). 

• In both these cases separation is forced at a sharp edge. There 
are several “boat tail” test cases (mainly axisymmetric) with sep- 
aration from a smooth surface, generally with small separation 
regions. 
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WHY MODEL TURBULENCE? 


. . .because, of course, we are too cheap to solve the exact equa- 
tions like the structures people do 

• "Turbulence modeling is the pacing item of CFD” - the late 
Jack Nielsen and many others 

• Stan Birch of Boeing said - “what limits you is computing 
power, not turbulence models” 

• Rightly or wrongly, industrial users stick with eddy-viscosity 
models 
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“BEST BUY” = WHAT I HAVE NOW 


• Correct - remembering that even testing of existing models can 
be limited by available/affordable computing 

• There is an urge to say "if my model does thin shear layers and 
my pet N-S case (even the backstep) it must be OK”. 

• Large eddy simulation just models the smallest (sub-grid-scale) 
eddies which do not carry much stress or momentum 

• But near a solid surface all eddies are small 
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EDDY VISCOSITY-1 


“Glushko jjses the TKE equation only to calculate the turbulent 
intensity q 2 and then makes the highly questionable assumption 
that the eddy viscosity vt is given by (q 2 ) 1 L j v? — constant" 
PB et al. JFM 28, 593 (1967) 

. . .(Glushko was implementing Prandtl-Wieghardt (1945)) 

• Eddy viscosity is easy to define (and measure). It is the ratio of 
a turbulence quantity to a mean-flow quantity, 

( Reynolds stress) /{Mean strain rate ) 

• Therefore it should not be treated as a pure turbulence quantity 
(like TKE say) 
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EDDY VISCOSITY-II 


• It may be better behaved / easier to correlate than Reynolds 
stress - but no guarantees 

• In self-similar flows with one velocity scale and one length scale 
(e.g. u T and 5) eddy viscosity must scale as 

vt = u r 8f(y/8 ) 

(and in slowly-changing flows near a solid surface f{y/8) = k y/8 
so v t = Ku r y) 

• In such flows, ( production / dissipation) = g{y/8 and near a 
solid surface g = 1 
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UNSTEADINESS - I 


• Reynolds averaging does not and cannot distinguish between 
unsteadiness and turbulence 

• Neither can any other sort of averaging unless the unsteadi- 
ness is periodic or otherwise structured AND occupies a different 
frequency, or wavelength, range from the turbulence. 

• Unsteadiness in a separated flow can very sensitive to boundary, 
or initial, conditions. 

• This makes me pessimistic about the reliability of (Reynolds-) 
averaged models in separated flow. . . 
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UNSTEADINESS - II 


. . .and of course we often want to predict the unsteadiness, to 
estimate structural buffeting or low-frequency noise 

• Flows which are very sensitive to boundary or initial conditions 
are ill-posed problems, and there’s an end of it! 

• Better-posed problems need 3D time-dependent simulations, 
DNS, LES, DES or whatever 

• DES, fed by a Reynolds-averaged model is likely to be reliable 
only if the boundary layers at separation are so thin that details 
of the turbulence don’t matter 
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MY TAKE ON REYNOLDS NUMBER 


• (1) Direct effects of viscosity on Reynolds stress in fully-turbulent 
flow are small for u T y/v < N where N is a function of space and 
time. N ~ 200 at I IT and KTH (at present!). 

• (2) However, the surface of the "viscous superlayer” has a di- 
mension much greater than 2 and its volume occupies a large 
fraction of the intermittent outer layer at low (i.e. laboratory) 
Re, leading to “direct viscous effects”. 

• (3) Indirect effects of viscosity can enter (e.g. flat-plate bound- 
ary layers) via d8/dx or du T /dx 

• At “high” Re, (1) and (2) are not a concern. The received 
wisdom is that neither is (3). . . 
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MY TWO CENTS ON THE LOG LAW-1. 

• The bare-bones derivation of the log law, due to Landau, is that 

if u = f(u T ,y,u ) 



. . . and we expect f -¥ 1 /k (say) as its argument tends to infinity. 

• Obviously a necessary, but maybe not sufficient, condition for 
the "if" to be true is t/<(5 

• Any other derivation is just the addition of bells and whistles to 
this dimensional argument. 
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MY TWO CENTS ON THE LOG LAW-II. . . 


• The power laws of Barenblatt and colleagues imply that / does 
not tend to a constant but continues to vary as an analytic func- 
tion of u r y/u 

• (i) I find it very difficult to believe that a second approximation 
to the log law will be another simple analytic function! 

• (ii) Also, direct effects of viscosity must be small for large 
enough u T yjv. 

• (iii) However (earlier slide) the effect of streamwise gradients of 
<5 and u T (indirect effects of viscosity) may extend to quite large 
Re (consider the uncertainty over the “exquisitely sensitive” wake 
parameter II). 
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MY TWO CENTS ON THE LOG LAW-3. . . 


• Panton has rightly pointed out that the regions of claimed va- 
lidity of the log law and the Barenblatt et al. power law are not 
the same (though the power law is supposed to overlap part of 
the log law). . . 

. . .but this does not change the above discussion of possible indi- 
rect Re effects 

• My advice to all concerned:- don’t do second-order math, on 
first-order physics! 
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CONCLUSIONS 


• The "High Reynolds Number” part of our title is hopefully just 
a warning to beware of low-Re effects in our test cases. . . 

. . .but the “Separated Flows” part of the title is a reminder that 
such flow combine unsteadiness and turbulence 

• ( Reynolds-)a veraged models can’t distinguish the two in general 

• If we want to predict the unsteadiness with any confidence we 
need ( x,y,z,t ) simulations 
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REWRITING HISTORY 


• Surgeon-General's warning. . . I added this after the meeting 
- which I unfortunately missed so I couldn’t give my views in 
discussion 

• I didn't say so explicitly, above, but my pessimistic views about 
the (un)reliability of Reynolds-averaged models imply the general 
conclusion of the Workshop. . . 

. . .that these models need very careful validation 

• A model user needs validation in the type of flow that he/she 
wants to calculate (since a model with fixed coefficients is unlikely 
to give results in all flows). 

• If this implies using different coefficients in different flows, that s 
bad science but may be good engineering 
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Some physical aspects of separated flows 
(separation on flat-plate-like objects) 



“In the spirit of a speculative research proposal” 

Preliminary remarks on: 
o LB methods 

o Physics of turbulence modeling 
o Vortex dynamics 
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Drag Coefficient over multiple bodies (PowerFLOW 3.1+3.2 v. Expt) 
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Cp 
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Cp Vs X for the 2D Bump 
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PuwerFLOW 


Three basic ingredients of the caienlatioBiS 

1. Discrete Boltzmann systems 


d t f(x; v; t) + V.Vf(x; v; t )= C 

Lattice gas reduces the phase space tremendously 
BC simpler, mesh refinement easy to implement 

2. Turbulence models 

k-e models: relaxation time in C varies locally to 
model “subgrid” dynamics 

3. Treatment near the waM 

First point is patched to a modified log-law: 

U s + = A log (y + /£(Vp)) + B 
^(Vp) = 1 + a (u.Vp) 
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Surface Cp (Alpha = 12) 



cm '•-OT-osjcoTj-mcor^ 
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x (dirnless) 



CA-0012 Alpha=U EXA Corp 




Exa Corp. Coefficient of Lift vs. Angle of Attack 

NACA 0012 Foil 
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Angle (deg) 
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A lowest order hypothesis 

TC = Y7Ta + (l-Y)TCs 

k = Y k a + (1- y) k s , 8 = Y S a + (1- Y) £ s 

Working forms of k and e are known for the attached 

case 

For instance, in the overlap region 

k a = u^/C,* 172 , 8 a = m T 3 /Ky 

What is the form in the separated region? 

k s oc l 2 (dU/dy) 2 , 8s = Vt(5U/dy) 2 

l is most likely obtained through the mean velocity 

Two things are needed: 

^Characterization of the separated region 
*The intermittency factor 
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Low-dimensional dynamics may be possible. 
What sort of low dimensional dynamics? 



uptfrttiA 



L = exp -(L-<L>), <L> = y 


Intermittency factor: 

A s p'5/x w , A ~ (y - yo) 
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(Only a very rough sUteMAtl) 




Tfc*n*v: 


Idealized Teniu: 


Iaviscid iji stability gives: 


' rol!e<U ub ’ structures 


Vortex sheet 


growth mtt 




Similar statements (with less experimental backing) can be made about 
other free shear flows (eg., wakes). Here wmt to eato the consequences 


The flow: 


ImviSCIO 

Idealized version: 


^ r'« /•* 


^ v*/ 


boumL*ry Uytr 

- flat ftate 

Varttx sheet 

- Mall 

Imaae VorttX 

sheet 


Invisdd instability, viscous effects and the stable coafjgBratkm: ? ? 
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The distribution of the Reynolds shear stress in (a) boundary layer 
and (b) channel (plane Poi seal lie) flows, evaluated from the numerically 
computed e i gens olut ions. The eigensolution data in (a) were from 
Jordinson (1970), and in (b> from Thomas (1953). 
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Figure 2: Plots of the turbulent shear stress r+ as a function of y+/Rl /2 («) across t 

“ “ Dd and (*> near its ^ Mgb Reynolds number experiments. The sources tor 

experimental daU are Antonia ct aL 1992 (channel, A. = 2S6, 916), Comte-Bellot 
1963 (c h a nne l, JL — 4324, 7309), Kim et aL 1987 (channel I>NS iL =- 1801 Laufer 

arj &.* s* n77 ' *">• i9 « w. rsis. 

1991 (channel DNS, A. = 125), Wei k Willmarth 1989 (dmnaeL R. = 71 £ 10201 
- d ZagaroU 1996 (pipe. A, = 851, 1430). The shear sSHS Ln otSnTby 
the numerical diSerentiation of the measured velocity ;,-oSe using Eq. (6) of section 
3. Zagarola's data for higher Reynolds numbers could ^ot be used because the mean 
velocity data have not been measured close enough to the wall. 
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1. PeAR REYNOLDS STRESS LOCATION IS LIKE A 

CRlCTICAL LATER 



2. THE CRITICAL LAYER 13 ESSENTIALLY A 

VORTEX SHEET 



153 




154 


Data from "peak values 


Up/Uo 



♦ Simpson’s boundary layer in adverse pressure gradient 
(all states of separation) 

Other symbols: Constant pressure boundary layers, pipes 
and channel flows 
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EXPERIMENTAL ISSUES RELATED TO TURBULENT HIGH 
REYNOLDS NUMBER SEPARATED FLOWS 

Roger Simpson 
Virginia Tech 

QUESTIONS 

What do we know? 

What are the problems with the data? 

What data are needed? 


NATURE OF SEPARATED FLOWS 

Nominal 2-D mean flows have highly 3-D instantaneous flows 

SOME COMMON MEAN 3-D FLOWS 

Juncture flows 
Cross-flow separation 
Vortical Separations 
Swept wing separations 

Some 3-D mean flows are complex and are difficult to calculate, 
even for lower Reynolds number cases. 
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DATA THAT WE NEED 


CASES WITH COMPLETE DOCUMENTATION 

• Detailed free-stream conditions 

• Surface pressure distributions 

• Skin friction magnitude and direction 

• Simultaneous 3-velocity component measurements 
with high sufficiently high spatial resolution 
for smallest turbulent scales; produces mean velocities, 
Reynolds stresses, and triple products 

• Temporal and spatial correlations for time and 
length scales (multipoint data) 

• Surface pressure fluctuations for vibrations/noise issues 

• Rate of strain tensor, vorticity tensor, dissipation rate 


157 


o¥€r airfoil 



158 


.... ' -A.** ' a: - 

tedkJitmit'' ~ f° 

h 


MATURE OF A SEPARATING TURBULENT BOUNDARY LAYER , 

*•* s 03 



THE SMALL HEAN BACKFLOW ID DENOTES INCIPIENT DETACHMENT; 
ITD DENOTES INTERMITTENT TRANSITORY DETACHMENT; D DENOTES 
DETACHMENT THE DASHED LINE DENOTES U - 0 LOCATIONS. 

FROM SIMPSON EI_AL (1981b) 


1. LARGE EDDIES GROW DURING DETACHMENT) I 6 (ft SpOMUftSC ftiut. 

2. LARGE EDDIES SUPPLY TURBULENCE ENERGY TO BACKFLOW AND CONTROL 
OUTER REGION ENTRAINMENT RATE 

3. LARGE EDDY BEHAVIOR SCALES ON MAXIMUM SHEAR STRESS. 

A. DIFFUSION AND DISSIPATION OF TURBULENCE ENER6Y IN BACKFLOW. 

5. SMALL -uv IN BACKFLOW; COLES 'LAW-OF-THE-WAU/ DOES NOT APPLY. 

5nwps#R> Tvs z d* £ &**** t ^ tgukalfa t 

■//<?«*! A6Akb«J**l>k 287 ; 

Simpson (.mb') T>nj.Aetosf**' St * , 
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Aspects of turbulent boundary-layer separation 
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Fig, 22 Detaching flow mean velocity profiles (Slap son etal. , 1981a) 
from lastr and hot-wire anenMeter data. Note displaced ordinates. 
Sol'd line denotes Re e - 47000, H * 5.38 data of Hastings and Moreton 
(1982) 
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MeaM 
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U/U„= 0.3( Y/N - LN(Y/N) - 1) - 1 


Fig. 23 Normalized backflow wean velocity 
profiles. Simpson et al . (1981a). L8V: 0 • 
3.53 m. ; P* 3.68 mTo. 3.97 V* 4.34 «. 
Hastings and Mare ton (1962); LDV O. 

Uestphal pul sed-wirt • . Solid line, equation 
(IV. 11). Dashed line, equation (IV. 12). Fra 
Slnpson (1983). 
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AtKOSPAfcl! ANtfBfc* 



Virginia Polytechnic Institute and State University 


VOrtical separation o* Vump 



Perspective view of leeside (above) 

Surface Oilflow over Large Bump #3 
U = 90.2fps ,27.5mps), Reg = 6000, H/5 = 2, 5 = 1.54 in. = 39mm 

Top view (below) 
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Vector Length: Grid/M agnitude = 1.2 


Note: 

• The logarithmic 
scale for Y/H 

• The nearest wall 
dau at about 
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Figure 3: VW mean secondary flow vectors for 650 individual LDV data points at 
11.14” downstream of the center of Bump#3 with 5/H=l/2 (X/H=3.63V . 

Ohc xtreamcvt'se vorftyt we or 


£ 
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Figure 4: Typical kKO model computational solution for 
Bump#3 with 5/H=l/2 at X/H=6.5. 
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OIL FLOW PATTERN (COURTESY OF ARA, BEDFORD) 



INTERPRETATION OF SKIN-FRICTION LINES ON TOP 
SURFACE OF PORT WING 


(c) TURBULENT FLOW ON SWEPT WING OF FIGHTER 
AT HIGH LIFT 


| FLOW 



(d) TURBULENT FLOW ON INBOARD TRAILING EDGE 
REGION OF LIFTING HIGH ASPECT RATIO SWEPT 
WING OF TRANSPORT AIRCRAFT AT HIGH SUB 
SONIC SPEED (COURTESY OF B. ELSENAAR. NLR) 


Fig. 14 Continued. 

AGAflDoin^iC 
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SOME DIFFICULTIES OF TURBULENCE MODELING FOR 
THREE-DIMENSIONAL FLOWS 
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Filiations for Modeling 

x- direction momentum equation: 

DU „dU „dU-„ T dU -ldP _d 2 U duv 


Dt 


= U- 


dx 


W— = 


f v- 


dy dz p dx dy dy 


z- direction momentum equation: 

DW „dW dW -IdP _d 2 W dvw 


Dt 


= U- 


H V — r rr — -v 

dx dy dz p dz dy ' dy 


W 


Continuity Eq.; 
dU dV dW 

dx dy dz 


= 0 


Transport of - uv Stress: Jf 
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Z)r 3y P 3*J 3y 

Transport of - vw Stress 
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GL, LRR, SSG, GL, LRR, 

FLT1, y*>15 y + >15 







CONCLUSIONS 


1) Low uncertainty data set above y + =3 at 30 stations are used to 
obtain the transport rate budgets at 4 stations. 

2) The pressure-diffusion term largely affects the transport rate 
budget of the v 2 , UV and VW stresses at all stations near wall. 

3) At Station 5, Separation and Vortex-core stations the distributions 
of the transport budget terms decrease progressively more steeply 
to approximately zero values than the 2-D flow station at lower y + 
locations. 

4) Anisotropic dissipation rate model of Hallback et. al. shows that 
the dissipation rate is increasingly anisotropic with increased three- 
dimensionality. 

For 3-D flow stations dissipation rate distributions are 
approximately isotropic above y + 100. However for 2 -D flow it 
becomes isotropic above y + 350. 

5) Comparison of the present data to the DNS results show that the 
present data are in better agreement with the DNS solutions using 
the anisotropic dissipation rate approximation. 
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CONSTANT TEMPERATURE HOT-WIRE ANEMOMEETRY 

FOR 3-D FLOWS 


Advantages: 

• Continuous time series data from multiple sensors 

• Frequency response: 30KHz flat response; 
carefully matched phasing flat up to 50KHz; 
calibrated matched phasing up to 150KHz. 

Disadvantages: 

• +/- 40 degrees total included angle for 4 sensor probe for mean velocities 

• +/- 20 degrees for turbulence for triple products 

• 0.5mm smallest measurement volume 

• 1mm closest near- wall measurement 
• Near-wall flow interference, even with single wires 
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LASER DOPPLER VELOCIMETRY 


Advantages: 

• Linear Doppler frequency - velocity relationship 

• Miniature 3-velocity-compponent probes with low 
uncertainties of each velocity component 

• Fine spatial resolution measurements - 
30im spherical measurement volume within 30im of wall 

• Rate of strain, vorticity, dissipation measurements possible 

Disadvantages: 

• Flow seeding required 

• Low data rate - 10 3 to 10 4 coincident signals per sec. 

• Single point data 
• Setup and data acquisition time 
• Probe hardware restricted to model interior 
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TABLE 3. 20:1 odds +2 C uncertainties of 
means velocities, Reynolds’ stresses and 
triple products. 



U Ju x 
V/u T 
\N/U X 







uv/u. 



0.148 

TFv/u if 

0.082 

0.033 

u z w/u? 

0.144 

0.097 

v z w/u? 

0.051 

0.062 

TnF/u? 

0.055 

0.04 


0.165 

0.07 


0.053 

0.025 

UVW/l if 

0.05 

0.037 

tF/t / 3 

0.254 

0.019 

aa 

0.043 


wPlu? 

0.424 
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Receiving optics for AU/Ai, AU/Ay, AU/Az 



Figure 1. Schematic of AU and AV incident and receiving optics for ROSVOR 
LDV. AW incident and receiving optics same as arrangement for AU, 
but lying in the YZ plane. 
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PARTICLE IMAGE VELOCIMETRY 


Advantages: 

• Global measurements 

• Faster data acquisition time for one plane 

Disadvantages: 

• Flow seeding required 

• Setup and data processing time 

• Only planar data with lower uncertainties 
(out of plane data with much higher uncertainties) 

• Higher uncertainties for u\ v’ than hot-wire and LDV; 
much higher triple product uncertainties 

• Multiple fields of view required to resolve various scales 
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